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Abstract 

We investigate propagators in Lorentz (or Landau) gauge by Monte Carlo 
simulations. In order to be able to compare with perturbative calculations we 
use large /? values. There the breaking of the Z{2) symmetry turns out to be 
important for all of the four lattice directions. Therefore we make sure that the 
analysis is performed in the correct state. We discus implications of the gauge 
fixing mechanism and point out the form of the weak-coupling behavior to be 
expected in the presence of zero-momentum modes. Our numerical result is 
that the gluon propagator in the weak-coupling limit is strongly affected by 
zero-momentum modes. This is corroborated in detail by our quantitative 
comparison with analytical calculations. 



1. Introduction 



Fundamental results on QCD have been obtained in continuum perturbation theory 
as well as nonperturbatively on the lattice. It appears important to connect the 
knowledge from these two approaches more closely. Progress within this respect is 
expected from the calculations of typical quantities of perturbation theory by Monte 
Carlo simulations in lattice gauge theory. The gluon propagator is a fundamental 
quantity for such an investigation. 

In the seminal paper by Mandula and Ogilvie ^ the zero-momentum gluon prop- 
agator r(r) = r(r;p = 0) has been numerically calculated in the Lorentz (or 
Landau) gauge. The analysis of the r-dependence of this propagator drove the 
authors of to the conclusion that r(r) is consistent with the propagation of 
a massive gluon. Therefore, the effect of dynamical mass generation was declared 
to be observed. Since that there was quite a number of papers devoted to the 
calculation and analysis of gauge-variant zero-momentum propagators at zero and 
nonzero temperatures following the same hne as paper (see, e.g. [2-8]) However, 
the impact of zero-momentum modes on the propagators has not been analyzed. 

Recently it has been shown by one of us that zero-momentum modes may strongly 
affect some gauge- dependent correlators 0. In some cases, e.g., in the Coulomb 
phase in pure U{1) gauge theory, zero-momentum modes can mimic the effective 
masses defined in a standard way (see eq. ( |5.2| )). This is apparent from the data 



of |]T0|. Because so far nothing is known about the magnitude of this effect in 
nonabelian gauge theory, it appears highly desirable to perform simulations which 
clarify this issue. In order to be able to compare quantitatively with perturbative 
calculations and to draw firm conclusions it is reasonable to begin this study in the 
truly perturbative region, i.e. at large values of p. 

In the present paper we perform Monte Carlo simulations in pure SU{2) gauge 
theory in the four-dimensional volume V4 = L1L2L3L4 with periodic boundary 
conditions. We use the Wilson action 

^ = /?E(l-^Tr[/p) , (1.1) 



p 



and the Lorentz gauge as described in Section 2. We discuss the working of the gauge 
fixing mechanism in detail and present analytical forms of the propagators appro- 
priate at weak-coupling in the presence of zero-momentum modes. We demonstrate 
the effects of the broken Z{2) symmetry and show that in the appropriate Z{2) 
state there is agreement with perturbative calculations. In particular, we find that 
there is a sizable effect caused by zero-momentum modes which can be understood 
quantitatively. 

2. Gauge fixing procedure 

In SU (2) lattice gauge theory the Lorentz gauge is fixed by maximizing 




by appropriate gauge transformations U^^x yx+eJJfj.xV^ , which implies the local 
condition 

E^Mf^M- = (2-2) 

(where (9^/(x) = f{x) — f{x — e^)). The maximized functional ( |2.1| ) may be cast 
into the form 

^max = ^ETr (M.l-J) (2.3) 

where 

= E(K.+e,t/;.. + K-e,f/i.-eJ ■ (2-4) 

From the fact that ( p.3| ) is maximal it follows that the maximizing gauge transfor- 
mations are given by 

V, = MJ\\M4 (2.5) 

where ||Mj,|p = |Tr(Mj;Mj). Obviously (p^), is a system of equations which de- 
termines all V^s given the f/'s of the configuration. (It should not to be mixed up 
with the relation without V^s on the r.h.s. and with the U's being the ones of the 



3 



particular iteration step, which is used iteratively site by site in numerical gauge 
fixing) . 

The gauge-fixing procedure is based on the equality (P) / = (P^) of the expec- 
tation with gauge fixing for the observable P to that without gauge fixing for the 
transformed observable . We note that this actually amounts to switching to 
the gauge-invariant observable Pcfr which is related to the gauge-variant one P in 
the gauge under consideration, i.e. that Peg = P^ . The respective general relation 
(P) / = (Pcff) with Pcff = /y JFP/ /y JF , where JF is a general gauge fixing function. 



has been pointed out some time ago |11 



It should be realized that the mechanism leading from a gauge-variant observable 
P to an effectively gauge-invariant one Pefj is simple. The gauge-fixing function 
provides strings of gauge fields which can combine with the given gauge-variant 
operators to gauge-invariant objects (and, of course, only such objects contribute to 
the integrals). For example, given a gauge field at a certain link, a string may be 
provided such that one gets a loop which contains this link. In practice there can 
be many suitable strings such that one arrives at a sum of gauge-invariant terms. 

In the present case the gauge-invariant observables can be constructed. To see this 
one has to realize that after maximization the original U variables due to the gauge 
transformation are decorated by V factors. On the other hand, these V factors are 
given here by ( p. 51 ), which therefore can be used to eliminate the V^s successively 
from P^. In detail, one replaces a by ( |2.5| ) which introduces next V^'s from the 
r.h.s of (|2.5|) , then replaces these y's again by (|2.5| ), and so on. This procedure 
terminates when all V^s are eliminated and one has arrived at an expression which 
contains only closed loops of gauge fields U . 

In the gauge-invariant observable thus obtained there are loops which include the 
gauge-variant parts of the gauge-variant observable P according to the mechanism 



explained above. Furthermore, noting that the factors ||Ma;|| in (|2.5| ) also contain 
\^'s, it becomes obvious that the terms of the effectively invariant observable are 
multiplied by weights which stem from these \\Mx\\ factors and which only depend 
on gauge- field loops (the locations of which are not restricted). 
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This construction allows insight into properties to be expected. First of all it is 
obvious that very large loops and, in particular, also loops of Polyakov type occur. 
Thus one can expect that the effective gauge-invariant observables resulting from 
gauge- variant ones become sensitive to the breaking of the Z{2) symmetry. The 
gauge-variant propagators investigated here indeed turn out to be strongly affected 
by the occurring Z{2) states as will be discussed in Section 4. 

In the case of the propagators defined in Section 3 P contains gauge-variant 
terms of forms Ti^U^^Uiyy), 1i{U^JJly) and their hermitian conjugates. By the 
maximizing gauge transformation then, for example, the term Ti{U^xUyy) gets the 
form liiVx+eJJ iixV^Vy+eJJ fiyV^) . The elimination procedure of V"'s described above 
forms loops of gauge fields incorporating the U factors of these terms. In the non- 
abelian case because of the product form of these terms one necessarily has to 
connect x and y in order to eliminate all y s. In other words, one has to form loops 
connecting x and y while separate loops through x and y are not possible. This 
has the consequence that one actually does not have correlations between separate 
operators, which implies that the usual transfer-matrix considerations of spectral 
properties do not apply. 

We have performed Monte Carlo simulations on lattices of sizes 

43 

X 8 , 8^ X 16 and 

16^ X 32. After studies at various (3 the investigations to be described in Section 5 
have been at large (3 (mainly ai (3 = 10, some at /5 = 99). Similarly, after studies of 
transition rates between Z{2) symmetry states with heat bath as well as Metropolis 
algorithms, for the main simulations heat bath update has been used. Measurements 
have been usually separated by 100 sweeps, on the small lattice also by 10 sweeps. 
Gauge fixing has been done iteratively site by site using the local solution ( p. 41) 
(with no K's on the r.h.s. and with the ?7's being the ones of the particular iteration 



step) and applying stochastic overrelaxation . This iterative procedure has been 
stopped if ^E. II i:,d-0,{x)f < 10-8. 

Statistical errors of observables have been determined by the jackknife method. To 
control the quality of gauge fixing 10 gauge copies have been generated by random 
gauge transformations for the configurations used in the measurements and the data 
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separately analyzed for the best (largest F), the first and the worst (smallest F) 
copy. The differences have been found to be small as compared to the statistical 
errors of the observables of interest, which agrees with similar observations in other 
investigations |], |13|. Furthermore, it has also been checked that T^^r) is constant 
within errors as it should be. 

3. Forms of propagators 

Let us define traceless fields 0^{x) 

0,ix) = ^(t/^. - Up , U,, = e^^^- , (3.1) 
and consider propagators 



T,{p,r) = 5]Tr(a^(p,t + r)a^(-p,t)) (3.2) 
where (with V3 = L1L2L3) 



d,{p,r) = -j:e''^-'0,{x,T) . (3.3) 

We choose p = (0,0,^3) (where ^3 = with integer p) so that the transverse 
components are Fi and F2 (in the following we use Fy = |(Fi + F2) and analogous 
expressions for related quantities). 

Perturbative results serve usually as a reference form in the numerical (nonper- 
turbative) study. However, a perturbative approach in nonabelian gauge theory 
enclosed in a periodic box faces a problem |]T1|, |r3|. This is the zero mode problem 
which makes a normal perturbative expansion invalid because of the divergent gaus- 
sian integrals. These zero modes (i.e. zero-eigenvalue eigenstates of the quadratic 
form) are associated with zero momentum, so this problem is usually referred to as 
zero-momentum mode problem. 
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Ignoring the divergency of the gaussian integrals (which is usually the case) one 
arrives at the following expression for the free gluon propagator: 



3(^2 g-«P4-r 

r.(p, = ^ E 4sin^f + 4sin^f " ^''^^ 
Perturbatively the gluon is massless, and for p = this correlator is ill-defined. 

A naive implementation of the idea of dynamical mass generation suggests to 
generalize eq. (|3.4|) as 



^ 21^ 5 4 sin2 f + 4 sin2 f + A sinh^ ^ ' ^^'^^ 

irig being a 'gluon mass' (defined such that E = rrig for p = ). It is usual 
practice to interprete all deviations from eq. ( |3.4| ) as an effect of nonzero shown 
in eq. ( p.5| ) (see, e.g. In this work we would like to discuss an alternative to 

eq.'s O and (|^) for p = 0. 

There is a degeneracy of the x-independent (zero-momentum) solutions U^^ of 
the classical equations of motion due to the toroidal structure of the periodic lattice 
TB|. An example of such a solution with zero action (toron) is Uf^ = e*'^''^ , where 
T is one of the generators of the gauge group and 0^ are four numbers. The 
perturbative expansion deals with fluctuations about the classical solutions (e.g. 
torons), and the shift A"^^ ^ 0^ + A"^^ (for small fields) produces the appearance 
of a non-negative constant term ~ (0^) in the correlator r^(r). 



To get our alternative description to eq.'s ( p^) and ( |3.5|) we split off the zero- 
momentum parts 



^4 



X 



of the operators 0^(x) so that they decompose as 

d,{0,T) = C^ + 5d^{T) . (3.7) 
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Then eq. ( p.2| ) for p = gets the form 

r,(r) = Tt{cI) + R,{t); 



(3.8) 



R,ir) = ^Y.'^T{6d,{t + r)6d,{t)) . (3.9) 



Evidently, 



f Er.M = Tr(C^). (3.10) 



The first term on the r.h.s. of ( p.8|) , being defined in terms of the original fields, 
can be readily determined in the Monte Carlo simulations. The term -R^(t) in ( p. 81 ) 
describes correlations between fields 60^{t) with zero- momentum parts subtracted. 
Making the change of variables A^^ — > 0/^ + ^x^j and using a collective-coordinate 
method^ one obtains in gaussian approximation 

Mr) = + |£ E ^ . (3.11) 

which will be our reference form in the numerical study of the propagators r^(r). 



4. Broken Z(2) symmetry 

To monitor the breaking of the Z{2) symmetry we have measured Polyakov loops 

P. = ^T.T^^UU,.. (4.1) 

Because of the two possibilities for each direction (i.e., > or < 0) there are 

2^ = 16 states : (+ + ++), (+ + H — ), . . . , ( ). Fig. 1 shows the time history 

of the P^, Tt and r4 obtained on the 4^ x 8 lattice. It is seen that the are, in fact, 
strongly affected by the indicated states and, in particular, that also the space-like 
^which works here, of course, only in lowest approximation 
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directions are important within this respect. The close relation between and r4 
is also obvious. 

On larger lattices the same features are observed (though getting accurate data 
for illustrations by time histories gets less easy). The phenomenon occurs for all f3 
in the deconfinement region. The transition rate between the states is smaller on 
larger lattices and for larger (3. It depends, of course, on the algorithm which is 
used. For example, with heat bath updates there are considerably more transitions 
than with Metropolis. 

From the numerical analysis we find in more detail that the observables, depending 
on their symmetry properties, take different values in the states. Huge differences as 
well as smaller differences of the values taken by the observables in different states 
are observed. This is seen in Fig. 2 for the example of Fy at finite p on the 16'^ x 32 
lattice (for clarity of the figure the statistical errors are not shown; they range from 
less than 2 % for the (+ + ++)-state to about 20 % for states with the largest 
magnitudes of the propagators). 

It should be obvious that to obtain sensible results measurements of observables 
are to be done separately for each state. In particular, for comparison with usual 
perturbation theory only the (+ + ++)-state is appropriate. To improve statistics, 
instead of only selecting the (+ + ++)-data from the time history, the other config- 
urations (solely) for the purpose of determining observables may be transformed to 
(+ + ++). We have checked that this leads (faster) to the same results. The Z(2) 
transformation used, which for a particular direction leads from — to + , is simply a 
factor —1 applied to the field variables of this direction (for the time direction this 
has already been used in Ref. 0). 

5. Zero momentum modes 

We now concentrate on the (+ + ++)-state at large /5, i.e. (3 = 10 and /3 = 99. 
To test the accuracy of the lowest order calculations we have compared numerical 
results for ps with (|3^) . The deviations are generally smaller than about 10 % 
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and become smaller for larger (3 and for larger p^. They are certainly very small as 
compared to the effects to be discussed below. Thus it appears that for the present 
study we are sufficiently deep in the perturbative region. 



For p = the comparison with ( 3.11 ) requires the knowledge of (C^) which we 
also determine in the Monte Carlo simulations. In addition, we obtain important 
properties of by comparing the time histories of C^, 50^{t) and of 

X,(p) = Re-^^e^^^-^0,(x) (5.1) 

with p = (0,0, l^p) and p ^ 0, i.e. of an example of a nonzero-momentum part of 
the fields. 

Fig. 3 shows time histories of components C^, 50^ and (where = 
Ealc^^/S)*^;^ etc.). For 501, ^^(1) ^^(2) uniform Monte Carlo noise is seen. 
Its magnitude is larger if the three-momentum involved in the particular quantity 
gets smaller. The behavior of in Fig. 3 presented for r = 0, is the same for all 
r. For in addition to noise with magnitude comparable to that of 5(9^, surpris- 
ingly large variations are observed. These variations exhibit different patterns for 
different /i, while for different a we find the same pattern. 



A very interesting aspect of these results is that in the decomposition (|3.7| ) of the 
fields the parts and 50^ obviously behave quite differently. The observed 
large variations of appear to be a characteristic consequence of zero-momentum 
modes. They seem to be related to changes between different constant classical 
solutions as one may consider in the weak-coupling limit. 

Because of the considerable size of the large variations of it is to be expected 
that also (C^) gets large as compared to other quantities. This will be seen to 
be indeed the case below. A further consequence of these variations is that error 
analyses of observables involving give unusually large errors. This, in particular 
holds for the propagator (|3.8|) . In order to show that the unusually large errors 
are entirely due to the part Tr(C^) we also have measured R^ij) separately in the 
simulations. An example of the respective results is shown in Fig. 4. It is seen that 
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while Tt{0, t) exhibits very large errors, for Rt{t) one indeed obtains errors of usual 
size. 

From Fig. 4 it is obvious that our numerical results for the transverse propagator 



with p = agree reasonably well with that of our lowest-order calculation ( |3.11 
The average over all times of the propagator is precisely Tr(C^) as predicted by 
(13.101) . Analogous observations as on the 8^ x 16 lattice for /3 = 10 have been made 



on the 4^ X 8 lattice for /3 = 10 and on the 16^ x 32 lattice for /3 = 10 and {3 = 99. 
Thus our numerical results confirm the description we have given. 

To quantify the importance of zero-momentum modes we compare Tr(C^) with 
i?^(0) (Table 1) and also with r^(p, 0) with p ^ (Table 2). To allow for the 
comparison of different lattice sizes and different (3 the values in Tables 1 and 2 are 
multiplied by V^/g"^ (which obviously cancels the extra factor implicit in our defini- 
tions as compared to usual continuum expressions). While all the other quantities 
then are roughly of the same order of magnitude, the values of Tr(C^) turn out to 
be much larger. From Table 1 it is seen that for increasing f3 this feature gets even 
more pronounced. The same holds for increasing lattice size. Not only ( gets larger 
but there is also an increase of the ratios (/j. 

In literature investigations of the gluon mass have been based on determinations 
of effective masses m(r) from 

cosh(m(r)(r + l-f)) ^ Tt{0,t+1) 
cosh(m(r)(r- ^)) 1^(0, r) 

For the present data the respective results are depicted in Figure 5. Because m{T)L^ 
shows little dependence on lattice size, these results may also be considered from the 
point of view of finite temperatures where screening masses are determined ||^, 
We have shown that instead of referring to effective masses the data can perfectly 
be described by eq. ( ^.llj ) relying on zero-momentum modes. 
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6. Discussion and conclusions 



To summarize, the dependence of the gauge-variant (zero momentum) propagators 
r^(r) on zero-momentum modes can be very strong. The mechanism of this influ- 
ence can be viewed as follows. The chosen gauge fixing does not prohibit the fluctu- 
ations about some nonzero constant values 0^ 7^ , and the shift A";.^ 0^ + ^xfi 
(for small fields) produces the appearence of a non-negative constant term in the 
correlator r^(r) Const + R^ij) . In this paper we demonstrated the occurrence 
of a large constant term explicitly in the weak coupling (perturbative) limit. 

Considering time histories in our simulations we have found that the zero- 
momentum term shows unusually large fluctuations. These fluctuations appear to 
be a characteristic signal for zero-momentum modes. 

The question of the role of the zero-momentum modes at smaller values of (5 
(in the physical region) needs further study. Indeed, one expects dynamical mass 
generation in the physical (nonperturbative) regime. However, zero-momentum 
modes which appear in the periodic volume may (at least, partially) mimic it, as it 
happens in the pure gauge U{1) theory. The latter can be seen from the data in 
|10| as has been discussed in To disentangle these two effects, i.e. dynamically 



generated mass and zero-momentum mode contribution, a comparative study with 
different boundary conditions could be of use. 

On the finite lattice at large (3 values one is beyond the finite-volume transitions, 
and thus the propagators become sensitive to the broken Z{2) symmetry states of 
the deconfinement region usually monitored by Polyakov loops. We find here that 
the states signalled by the space-like Polyakov loops are important, too. 
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Table 1 

C = {V^lg^)TT{{CTf) and 7 = {V^l g^){{Rm - Rt{L,/2)) 



lattice 




c 


7 


43 X 8 


10 


6.30 (16) 


1.610 (5) 


8^ X 16 


10 


15.7 (7) 


3.229 (17) 


16-'^ X 32 


10 


57 (10) 


6.03 (22) 


153 X 32 


99 


375 (86) 


4.98 (19) 



Table 2 



{Vs/g')rT{p, 0) with p = (0, 0, gp) 



lattice 


/5 


p=l 


p = 2 


p = 4 


43 X 8 


10 


0.460 (2) 


0.275 (1) 




8^ X 16 


10 


0.992 (4) 


0.457 (1) 


0.277 (1) 


16^ X 32 


10 


2.100 (20) 


0.981 (4) 


0.457 (1) 


16^ X 32 


99 


2.002 (23) 


0.926 (5) 


0.436 (2) 
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Figure captions 



Fig. 1. Time history of Polyakov loops and propagators (at r = 0) 

for P = 10 and 4^ x 8 lattice. 
Fig. 2. Transverse propagators in various states for p = (0, 0, |^) , 

(3 = 10 and 16^ x 32 lattice (curve is (^). 
Fig. 3. Time histories of C^, SO"^, X^{1) and X^(2) 

for (3 = 10 and 16^ x 32 lattice. 
Fig. 4. rr(0, t) and Rt{t) compared with ( p.ll| ) (curves) and Tr(C^) 

(constant line in upper Figure) for (3 = 10 and 8^ x 16 lattice. 

Errors of Rt are smaller than symbols. 
Fig. 5. Effective masses from eq. (|5.2|) for p = 0, (3 = 10 and 

various lattice sizes (curve based on ( p.ll| )). 
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Figure 1: Time history of Polyakov loops and propagators at r = for /3 = 10 and 
4^ X 8 lattice. 
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Figure 2: Transverse propagators in various states for p = (0, 0, , /? = 10 and 
16^ X 32 lattice (curve is Q). 
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Figure 3: Time histories of C^, 501, ^m(I) ^^^(2) /3 = 10 and 16=^ x 32 
lattice. 
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Figure 5: Effective masses from eq. ( |5.2| ) for p = 0, /3 = 10 and various lattice sizes 
(curve based on ( |3.11D ). 
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